Inertia and chiral edge modes of a skyrmion magnetic bubble 
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The dynamics of a vortex in a thin-film ferromagnet resembles the motion of a charged massless 
particle in a uniform magnetic field. Similar dynamics is expected for other magnetic textures with a 
nonzero skyrmion number. However, recent numerical simulations revealed that skyrmion magnetic 
bubbles show significant deviations from this model. We show that a skyrmion bubble possesses 
inertia and derive its mass from the standard theory of a thin-film ferromagnet. Besides center-of- 
mass motion, other low energy modes are waves on the edge of the bubble traveling with different 
speeds in opposite directions. 



Dynamics of topological defects is a topic of long- 
standing interest in magnetism. The attention to it stems 
from rich basic physics as well as from its connection to 
technological applications Theory of magnetization 
dynamics in ferromagnets well below the critical temper- 
ature is based on the Landau-Lifshitz equation 2] for the 
unit vector of magnetization m(r) = M(r)/Af, 



rii — 7B X m + am x rh, 



(1) 



where 7 is the gyromagnetic ratio, a <^ 1 is a phe- 
nomenological damping constant [3(, and the effective 
magnetic field is a functional derivative of the free en- 
ergy, B(r) = — 5[//(5M(r). The latter includes local 
(e.g., exchange and anisotropy) as well as long-range 
(dipolar) interactions, thus making Eq. ([1]) a nonlinear 
and nonlocal partial differential equation with multiple 
length and time scales solvable in only a few simple cases. 
For example, translational motion of a rigid texture, 
m — m(r — R(i)), is fully parametrized by the texture's 
"center of mass" R. For steady motion, R(t) = Vt, 
the velocity can be obtained from Thiele's equation Q 
expressing the balance of gyrotropic, conservative, and 
viscous forces: 




GxR + F-i:>R = 0. 



(2) 



Here G is a gyrocoupling vector, F — —dU/dR is the net 
conservative force, and _D is a dissipation tensor. A rigid 
texture moves like a massless particle with electric charge 
in a magnetic field and an external potential through a 
viscous medium. If G 7^ 0, the "Lorentz force" greatly 
exceeds the viscous drag. We thus ignore dissipation. 

Although Eq. 1^ was derived for steady motion, Thiele 
anticipated that it could serve as a good first approxi- 
mation in more general situations. Indeed, his equation 
describes very well the dynamics of vortices in thin ferro- 
magnetic films 043 ■ t'^is case, the gyrocoupling vec- 
tor G = {0,0, G) is proportional to a topological invari- 
ant known as the skyrmion charge q — (1/ Air) J dxdym- 
{dxin X dym), the film thickness t, and the density of an- 
gular momentum M/7; to wit, G = AnqtM/"/. A vortex 
has q = ±1/2 and thus G ^ 0. In a parabolic potential 
weh, U{X, Y) = JC{X^ + Y^)/2, it moves in a circle at a 
frequency uj = K/Q. 



FIG. 1. A circular magnetic bubble with skyrmion number 
q = —1. Out-of-plane magnetization is blue for Mz > and 
red for < 0; in-plane magnetization on the domain wall is 
shown by arrows, (a) Equilibrium, (b) Elliptic deformation. 
The dashed line marks the equilibrium position of the domain 
wall, (c) The trajectory of the bubble's center observed by 
Moutafis et al. Dots mark positions evenly spaced in time, 
(d) A 5-cusped hypocycloid with its directing and generating 
circles. 



Similar behavior is expected for other topologically 
nontrivial textures, e.g., magnetic bubbles in thin films 
with magnetization normal to the plane of the film [lol - 
IT^ . A bubble is a circular domain with rUz < sur- 
rounded by a domain with ruz > 0, or vice versa, Fig .flfa) 

. Recently discovered skyrmion crystals llj 
ticularly those which are found in thin films [If 
are periodic arrays of magnetic bubbles with the same 
skyrmion charge q — ±1 Zang et al. [l^ modeled 

skyrmions in these structures as massless particles with 
Thiele's dynami cs (pi ). However, numerical simulations of 
Moutafis et al. [20| revealed that this model fails badly 
for an isolated skyrmion bubble. In a parabolic potential. 
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the trajectory of the bubble's center was not a circle, but 
"roughly a pentagon," Fig. [TJc). 

The main goal of this Letter is to derive the correct dy- 
namical model of a skyrmion magnetic bubble. We first 
introduce our model phenomenologically and then derive 
it from the standard theory of a thin- film fcrromagnet. 
This allows us to characterize not only the center-of-mass 
motion of the bubble but also the dynamics of its shape 
within the same framework. 

Phenomenology. The puzzling trajectory of the bub- 
ble's center is readily reproduced if we endow a skyrmion 
bubble with inertial mass M: 

-MR + gx 11-101 = 0, (3) 

The new equation of motion is of the second order in d/dt 
and has two circular modes with frequencies 

uj± = -g/2M ± ^{g/2My+]C/M. (4) 

A particle with zero initial velocity follows a hypocycloid. 
If uj^/lu^ = —4, the hypocycloid has 5 cusps and indeed 
resembles a pentagon. Fig. [Ijd) . We repeated the simu- 
lations of Moutafis et al. and found that the motion of 
the bubble is described with good accuracy by a superpo- 
sition of two underdamped modes with eigenfrequencies 
a;/27r = 0.97 GHz and -4.27 GHz [2l|. 

To understand the origin of inertia, we shift attention 
from the center of the bubble, where nothing is happen- 
ing, to its boundary, a domain wall defined as a line y{x) 
where Mz{x,y) = 0. A nearly circular domain wall is 
conveniently parametrized in polar coordinates (r, <j)) : 

r(^)=f + ^r™e""*. (5) 

m 

The Fourier amplitudes r„i = r* ,„ describe waves with 
wavenumbers k = m/f traveling along the circular edge: 
To is the breathing mode, ri = {X — iY)/2 encodes 
the location of the center of mass, r2 parametrizes el- 
liptical deformations. Fig. [Ijb), and so on. On the do- 
main wall, magnetization lies in the plane of the film, 
m = (cos "0, sin '0,0). For a circular wall in equilibrium, 
m points along the direction of the wall, ip = (f) ± 7r/2, 
Fig. [Tfa). More generally, 

1^(0) = 0±7r/2 + ^0„e™^ (6) 

m 

The fields r(0) and ■0(0) are coupled to each other, and 
so are their harmonics and ipm. Integrating out ■0i 
generates kinetic energy for the center of mass. 

Theory. We derive the dynamics of transverse fiuctua- 
tions of a Bloch domain wall, first for a straight wall and 
then for a circular one. To this end, we employ a method 
of collective coordinatesgeneralizing Thiele's approach 
beyond steady motion |22l |. The Lagrangian formalism 
allows us to easily integrate out the hidden degree of 
freedom — in-plane magnetization — in favor of the more 
evident transverse motion. 



An evolving magnetic texture m(r, t) can be 
parametrized by a (potentially infinite) set of collective 
coordinates $,{t) = {S,i[t) t S,2{t) , ■ ■ ■}■ Their equations of 
motion are similar to Thiele's equation ([2]): 

G.,0 + ^-A,6 = 0, (7) 

with generalized forces Fi = —dU/d^i, gyrotropic co- 
efficients Gij = ~Gji, and viscosity coefficients Dij = 
Dji. Eq. ([7]) can be obtained from a Lagrangian L = 
A • I — U{^), where A(|) is a gau ge p otential with cur- 
vature dAj/d^i — dAi/d^j = Gij [23] ■ The gauge term 
contributes to the action S = J Ldt a, time-independent 
piece J A ■ d^ known as Berry's geometric phase. 

Straight domain wall. We first consider the dynamics 
of a domain wall stretched along the x-axis, y{x,t) « 0, 
'il'{x,t) sa 0, from x = to x = £. Its Lagrangian, 

L[y,V]= / dxgytP-U[y,tP], (8) 
Jo 

contains a gauge term with gyrotropic coupling g = 
2tM/7 [13,12^. The resulting equations of motion are 
—gip — 6U /dy = 0, gy — SU /St/j = in the absence of dis- 
sipation. The in-planc magnetization is aligned with the 
wall in equilibrium, ip = y' = dy/dx; the cost of small 
misalignments is quadratic in ^0 — y', so 

L[y,^] = / dx [gyi; - «;(^ - y'f/2] ~ U[y], (9) 
Jo 

with the stiffness k to be determined. The field ip can 
be integrated out with the aid of its equation of motion, 
gy — k{iP — y') = 0, to obtain a Lagrangian for transverse 
displacements, 

L[y] = dx {py^/2 + gyy') - U\yl (10) 
Jo 

where p = g^ j n is the Doring mass density 0, H^. 

Potential energy of a domain wall can be split into 
local and long-range contributions. The local term is 
proportional to the length of the wall and tension a: 

Vi\y\ = \ G^fd^^Td^ ^Ui\0\+ \ dxay'^/2, (11) 

where Ui [0] = at is the energy of a straight wall. Thus 

L[y] = Jdx {pf/2 + gyy' - ay'^ 12) ~ U^i{y\. (12) 

Neglecting for the moment the nonlocal term [/„;[?/], we 
obtain a wave equation, py -\- 2gy' — ay" = 0, with waves 
traveling left and right at different speeds — cf. Eq. (|4]): 

w = p-^ [-gk ± ^/gW+Vl^^ . (13) 

In the limit a <C g^ / p = k, the slow wave has velocity 
vi ^ a/2g that is insensitive to inertia; in-plane mag- 
netization adiabatically aligns with the direction of the 
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wall. The fast mode with velocity V2 ~ —2g/p involves 
oscillations of in-plane magnetization out of phase with 
those of the direction of the wall. 

A minimal model of a thin-film ferromagnet with out- 
of-plane magnetization includes exchange coupling A and 
easy-axis anisotropy K strong enough to overcome the 
dipolar shape anisotropy: the "quality factor" Q = 
2K/iJ,()M^ must exceed 1. In this model [2l|, a domain 
wall has the width X/y/Q — 1, where A — y/2A/ hqAP is 
the exchange length, and tension cr — 2fioAPtXy/Q — 1. 
The coupling between the in-plane magnetization and the 
direction of the wall is k — <7/{Q — 1). By using the ma- 
terial parameters characteristic of FePt with thickness 
t = 32 nm [iil, we obtain vi = 370 m/s and V2 = -2100 
m/s in the local model. 

The nonlocal part of the potential energy comes from 
long-range dipolar interactions. A domain wall in a thin 
ferromagnetic film produces a stray magnetic field whose 
energy can be written as a double line integral along the 
domain wall l26l[27l. 



r 1 '^d f dri ■ dr2 



where ad — hqAI'^P /tt is "dipolar tension." This expres- 
sion diverges at both short and long length scales and 
thus requires both short and long-distance cutoffs (pro- 
vided by the film thickness and the wall length). The 
local potential energy (fTTj) may be absorbed into the 
nonlocal part (|14p at the expense of renormalizing the 
short-distance cutoff. Expanding Eq. ([14]) to the second 
order in y yields the following result: 

U,4y] = Uni[0] + i —<Jdk^ In {ka)\yk\\ (15) 

Jo 

where yk = r^'^ dxy{x)e-'''\ C/„z[0] ~ -ad\^{l/a) 
is the energy of a straight domain wall, a — 
(t/2)e'-^^^+°"/'^'' is a short-distance length scale, and 
C = 0.577. . . is the Euler constant. The wave stiffness 
CTrffc^ In (fca) is negative for k < 1/a, which means that 
a straight domain wall is unstable against small defor- 
mations. This is the fingering instability occurring in 
systems with long-range interactions |26:, 28]. It can be 
prevented by placing the domain wall in a strip of fi- 
nite width w. Repulsion from the edges, mediated by 
a stray magnetic field, increases the wave stiffness by a 
fc- independent term Aaa/w^ ■ The wave frequencies are 



uj = p 



-gk± y^g^k^ + aap {k^ In \ka\ + A/w^) 



(16) 

The frequency spectrum (Fig. [2]) has two significant 
changes from the local model (fT3t : a gap opens up; the 
band bottom is shifted to fco — e^^l'^a"^ . 

Magnetic bubble. The dynamics of a circular domain 
wall is derived along similar lines. For a wall of given 
shape ([S|), the in-plane magnetization tends to align it- 
self with the waU, Voq(</') = ± 7r/2 - f-^dr/d(j). The 
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FIG. 2. The spectrum of spin waves u)/2ti on a circular do- 
main wall: numerical simulation (circles), model calculations 
for a straight wall with k = m/f (line). 



V'-dependent terms in the Lagrangian density are thus 
gfil) - K{ip - il'cq)^/2. The La grangian of the Fourier 



(14) modes © and © is 



L[r,-0] = 27rr^ ( gfmtp. 



U\r] 



(17) 



Integrating out tp yields a Lagrangian for r containing 
both kinetic energy and a Berry-phase term: 

LM = (7rrp|rmP - Anmi g r*^rrn) - U[r]. (18) 



We thus arrive at a Lagrangian for the center-of-mass 
mode ri = {X - iY)/2: 

L{X,Y) ^ M{X^+Y^)/2+gXY-ICiX^+Y'^)/2. (19) 

It yields the anticipated inertial dynamics of a mag- 
netic bubble 13]). The gyrotropic constant, Q = 2ng — 
AntM/'f = 2.29 x lO^^^ J s/m^, depends only on the 
topology of the bubble (here the skyrmion number q = I) 
and on the area density of angular momentum tM/j; 
therefore, it can be taken at face value. The spring 
constant /C comes from magnetostatic repulsion between 
the domain wall and the edge of the disk; its calculated 
value, /Ccaic = 0.020 J/m^ for f = 37 nm, should also 
be reliable. The weakest link in our theory is the mass 
term Alcaic = Trrp = ■nfg'^/K = 5.4 x 10~^^ kg. It de- 
pends on the coupling k between in-plane magnetization 
and the direction of the wall and thus requires an ac- 
curate model of the domain wall. Alternatively, the ba- 
sic constant Q can be combined with the measured fre- 
quencies, a;sim/27r — 0.97 GHz and -4.27 GHz, to obtain 
/C,im = 0.018 J/m2 and TWsim = l-H x 10^22 
expected, we find a good match between the calculated 
and simulated values of the spring constant K, but the 
M values differ by a factor of 2. 

We have performed numerical simulations to measure 
the frequency spectrum of waves with higher azimuthal 
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numbers m. We used the same geometry and mate- 
rial parameters as Moutafis et ai: a FePt disk of ra- 
dius R = 80 nm, thickness t — 32 nm, magnetization 
M = 10^ A/m, exchange constant A = 10^^^ J/m, easy- 
axis anisotropy K — 1.3 x 10^ J/m"^, and gyromagnetic 
ratio 7 = 1.75 x 10^^ s A/kg. These give the quality fac- 
tor Q — 2.06 and exchange length A = 4.0 nm. We uti- 
lized micromagnetic simulator OOMMF ,29!] in the two- 
dimensional regime with a unit cell of 1.25 nm. The 
equilibrium radius of the bubble was f = 37 nm. The 
free motion of the mth harmonic of p{(j)) was fitted by 
a sum of two underdamped components. The extracted 
frequencies are shown in Fig. [2] The theoretical curve is 
the straight- wall spectrum ([T6l) with k = m/f. We used 
the mass density extracted from the simulation of the 
m = 1 mode, p — A^sim/Trf, and set the effective width 
w equal to the disk radius. The theory works quite well 
for the slow mode, less so for the fast one. 

Discussion. The widely used Thiele's equation ^ 
predicts that magnetic textures with a skyrmion num- 
ber q 7^ behave as massless particles moving in a 
uniform magnetic field and an external potential. Al- 



though this approach works very well for magnetic vor- 
tices {q — ±1/2), it fails for skyrmion magnetic bubbles 
{q = ±1). Here we have shown that a skyrmion bubble 
behaves as a massive object and have explained the ori- 
gin of its mass. A skyrmion bubble possesses additional 
modes, which are best viewed as transverse fluctuations 
of its edge. These waves are chiral, i.e., they propagate 
with different speeds in opposite directions. The non- 
reciprocal wave propagation is expected in other geome- 
tries, e.g., striped and labyrinthine domains. 

In a skyrmion crystal, the discrete modes of a bubble 
turn into excitation branches. The slow and fast m — 
1 modes give rise to the magnetophonon and cyclotron 
branches [30| ; the cyclotron frequency uj = G / is in the 
GHz range. The m = breathing mode has been seen in 
numerical simulations [31]. The breating mode and one 
of the m — 1 modes have been found in Cu20Se03 [3^ . 
Branches with higher m may also be detectable. 
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Appendix A: Center-of-mass motion 

The motion of the center of mass {X, Y) of a magnetic 
bubble is conveniently represented by a complex vari- 
able ri = {X — iY)/2. The general motion of a bubble 
with inertial mass M and a gyrotropic coefficient Q in a 
parabolic potential with stiffness /C is a superposition of 
two spiral motions: 



ri{t) ^ a„ 



(Al) 



m=±l 



For weak dissipation, Tm ^ (^m, 



UJ±l 




(A2) 

The top panel of Fig.[T]shows the best fit to Eq. (jAip with 
w+i/27r 0.97 GHz, T+i = 0.25 ns^^, w_i/27r = -4.27 
GHz, and r_i = 1.0 ns^^. 

A small systematic deviation between the best-fit line 
and the data is plotted in the bottom panel of Fig. [T] 
Much of it can be accounted for by a superposition of two 
spiral motions with frequencies w+i iwo, where a;o/27r = 
3.15 GHz is the eigenfrequency of the bubble's breathing 
mode. This is likely the effect of an anharmonic coupling 
between the two modes. 

We thus conclude that any unidentified modes con- 
tribute no more than 0.1 nm in amplitude, or about 1 




0.4 0.6 O.i 

time, ns 

FIG. 1. Top panel: The displacement of the bubble's center 
of mass {X, Y) as a function of time (points) and the best fit 
to the two- mode approximation (|A1|) (lines). Bottom panel: 
The difference between the data and the best fit (points) and 
the best fit to a superposition of two spiral motions with fre- 
quencies ij+i±LUo (lines). The vertical scales in the two panels 
differ by a factor of 30. 




(i+lj) 



FIG. 2. Tracing a domain wall (solid line). Circles denote the 
centers of lattice cells in the simulation. Squares are points 
at cell boundaries (dashed lines). Crosses are locations of 
the domain wall obtained by linear extrapolation. Open red 
symbols mean AIz > 0, filled blue symbols mean AIz < 0. 



per cent, to the center-of-mass motion. 



Appendix B: Locating the domain wall 

We define the domain wall as a line of points in 
the plane of the film where out-of-plane magnetization 
Mz{x,y) vanishes. Because our numerical simulations 
were done on a discrete lattice, we needed an algorithm 
to extract a continuous line from discrete data points. 
The location of the domain wall was determined in two 
steps: we first identified a discrete set of points with 
Mz — 0; we then fit their positions to a line, Eq. (5) in 
the main text. 

The first step is illustrated in Fig. [51 The magnetiza- 
tion at the boundary of a cell is determined by interpola- 
tion with the magnetization at the center of the current 
cell and its nearest neighbors as follows: 



M 



(M, J + M 



M 
M 



A 
A 

+ M 



M. 



)/4. 



If the wall intersects the cell then the sign of Mz at one 
or more points at its boundary differs from that at the 
center of the cell. At the crudest level, we could use the 
centers of intersected cells as a proxy for the location of 
the domain wall. To refine this result, we used a linear 
interpolation for Mz between the center of a cell and 
the points on its boundary to find the locations with 
Mz — 0. The refinement reduced the discretization noise 
by a factor of 70, which was particularly important for 
modes with higher azimuthal numbers, which had small 
amplitudes. 
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Appendix C: Local interactions 

Tension a and alignment stiffness k come from the local 
portion of the energy functional. Consider a ferromag- 
netic film of thickness t with two domains of magnetiza- 
tion, m = (0,0,1) and m = (0,0,-1), separated by a 
straight domain wall of length i along the x-axis. The 
magnetization field can be parametrized as 

m(j/) = (sin d cos -0, sin Q sin V', cos ff) 

with B = 9{y) and tp = const. In the state of lowest en- 
ergy, in-plane magnetization points along the wall, "0 = 
(Bloch wall). The energy contains three terms: exchange 



^exchange 

easy-axis anisotropy 



At J (fr{d,ni)idim) 



AU 



dy 



C/anisotropy = ~Kt / (fruTil = -KU / dy cos^ e, 



where A = yj2A/ ji^IvP is the exchange length and Q = 
2K/fiQM^ > 1 is the "quality factor." The minimized 
energy per unit length gives line tension 



(C5) 



When in-plane magnetization deviates from the direc- 
tion of the wall, it creates bulk magnetic charges with 
volume density — V-M = —dMy/dy. This induces an ad- 
ditional magnetic field Hy « —My and thus generates ad- 
ditional energy density fioHy /2 « sin^ sin^ ip. 
The domain-wall energy per unit length increases to 

U{'ip)/i = 2fioM^tX^jQ- 1 + sin2i/; - a + K0V2 (C6) 

for small angles 0. The strength of coupling between the 
azimuthal angle and the direction of the wall is 



1 d^U 



I rf02 



i/j=0 



Q-1 



(C7) 



and the energy of the magnetic field, which can be eval- 
uated as the Coulomb interaction of magnetic charges at 
the top and bottom surfaces of the film with area densi- 
ties tM^{y): 

f/dipoiar = £ / '^'^i ^^"^2 ^^(ri - r2)M,(ri)M,(r2) 



with the interaction kernel 
1 

r 



V{v) 



1 



Vr2 + i2 



V(r) 



For future convenience, we extend the definition of V(x) 
to negative values of the argument so that V(— a;) = N{x): 



1 



t2 



= lim 



1 



1 



^2 

(CI) 

The interaction kernel V{y) has a peak with a charac- 
teristic width of the order of the film thickness t. If 
magnetization varies slowly on that length scale, we may 
approximate V{y) k, 27r<5(r). Then the dipolar energy 
assumes a local form: 



^dipola 



t / d^r 



dy cos^ 9. 



(C2) 

Put another way, we make a local approximation for the 
magnetic field inside the film, i?z(r) ~ —Mz{r), whose 
energy density fioH'^/2 « fJ.oM'^/2. 

The energy cost of a domain wall is thus 



U ^£t 



dy 



AO' +{K - noM^/2) sin^ 



(C3) 



Minimization of this energy yields the domain- wall profile 

'y 



cos9{y) 



tanh ( ^ y/Q 



1 



(C4) 



Appendix D: Long-range interaction 

Equation (|C2p captures the local part of dipolar en- 
ergy, whose density is determined by the local value of 
magnetization. Inhomogeneities in magnetization pro- 
duce a stray magnetic field that gives rise to a nonlocal 
component. 



stray 



47r 



J d^n d\2 G(ri - rs) VM,(ri) • VM,(r2), 

(Dl) 

where r = (x^y) and V = (dx^dy) are two-dimensional 
vectors. The kernel is defined as a solution to the equa- 
tion V^G(r) = —V{r) and has the explicit form 

G(r) = ^^2 + t2 - r - t arsinh(t/r) = G(r). (D2) 

As we did for Y{x), wc extend the definition of G{x) to 
X < so that Gj{—x) = Gj{x). Note that 



'{x) + G'{x)/x = -V(x). 



(D3) 



For an infinitely sharp domain wall, Eq. (|Dip reduces 
to a double line integral, 

[/stray [r] - / dv, ■ dv^ G(ri - v^). (D4) 



The asymptotic form of the kernel (jD2[) is G(r) ~ —t^/2r 
ioT r t. Hence 



C/stray [r] 



with "dipolar tension" 



ad f dvi ■ dT2 



ri - r2 



2j-2 



(D5) 



(D6) 



7 



The simplified version of tlie stray-field interaction 
(|D5p is logarithmically divergent at short distances. One 
way to handle the divergence is to impose a short- 
distance cutoff, |ri — > b, where 6 is a length scale of 
the order of the film thickness t. The energy of a straight 
wall due to its stray field, computed with the exact ex- 
pression (jD4p . is 



C/stray ^ -^d [In i2£/t) + 1/2] , 

whereas the simplified version (jD5[) with a cutoff yields 

«7stray -od [In - 1] . (D7) 

The two expressions agree if we choose the cutoff to be 
6= (i/2)e-3/2. 

Local contributions to the energy of the domain wall 
ai can be absorbed into the stray-field energy (|D7|) at 
the expense of renormalizing the cutoff parameter. 



(D8) 



Appendix E: Energy of transverse fluctuations 

We compute the energy of a nearly straight domain 
wall, y{x) 0. The local part can be written as 



Jo ^ Jo 



dk 

2^ 



^k^\Vk\\ (El) 



where yu — dx y{x)e^^''^ is the spatial Fourier 

transform of y{x). 

The energy of the stray field (|D4[) can be expanded to 
0{y^) with the aid of the following results: 

dvi ■ dv2 = dxi dx2 [1 + y' {xi)y' {X2)] , 



The nonlocal portion of the energy of transverse fluctu- 
ations thus consists of two parts, 

C^stray[2/] = ^"^^ J dxi dx2 G{xi - X2)y' {xi)y' {X2) 

A^qM^ [ G'ix,~X2) , 2 

H / dxidx2—, ^{yi~y2) ■ 

n J 2[xi-X2) 

After an integration by parts and a Fourier transform, 
we obtain 



f°° dk 

[/stray [y] = ^ / :^^fc|yfc| 

Jo 2tt 



where 



111 Af^ f°° 

Ak = — / dx[l- cos {kx)] [G"{x) + G'ix) /x] 



With the aid of Eqs. ((CT|) and ((D3)) . we obtain 



[Koikt)+lii{kt/2) + C], 



where Kq{x) is a modified Bessel function and C = 
0.577 ... is the Euler constant. In the infrared limit, 
kt 0, 

Ak ^ <Jdk^ [\n(kt/2) + C - 1] . 
Upon adding the local term (jEl[) . we obtain 

dk 



U[y]^il ^cjd^\n{ka)\yu? 
Jo 2-K 

with the short-distance scale 

a ^ (t/2)e'^-i+'"/'^'^. 



(E2) 



(E3) 



